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Abstract 

We study the /?-deformed matrix models using the method of refined topological 
string theory. The refined holomorphic anomaly equation and boundary conditions 
near the singular divisors of the underlying geometry fix the refined amplitudes recur- 
sively. We provide exact test of the quantum integrality conjecture in the Nekrasov- 
Shatashvili limit. We check the higher genus exact formulae with perturbative matrix 
model calculations. 
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1 Introduction and Summary 





Recently there have been some interests in refined topological string theory, which origi- 
nate from the deformation in supersymmetric gauge theory [24j. It is expected that the 
Dijkgraaf-Vafa conjecture [8j, which relates matrix models with topological strings on cer- 
tain local Calabi-Yau manifolds, can be generalized to the refined case. Here the refinement 
will correspond to the /^-deformation of the matrix models. The topological expansion of 
ordinary matrix model free energy has been constructed from the spectral curves [21 [10] , 
and can be generalized to the /3-deformed case [H [7j . However, the topological recursion 
method is still too difficult for many practical calculations in the /3-deformed case, see e.g. 
[U [20]. Furthermore, the holomorphic anomaly equation can be derived from the topo- 
logical recursion in the undeformed case [9j. For the /3-deformed case such a derivation 
is not available at the moment in the literature. In this paper, we shall study the sim- 
ple example of /3-deformed cubic matrix model. Continuing in the direction of previous 
works \14:\ [161 [1] , we provide some higher genus formulae from the refined topological string 
method of holomorphic anomaly equation and the gap boundary conditions near singular 
divisors [3l[l8l[l5]. 



2 Perturbative calculations 

The partition function Z and free energy F of a Hermitian matrix model with polynomial 
potential 14/'(<I>) is defined by 

where gs is a perturbative expansion parameter, Aj (i = 1,2- •• ,iV) are the eigenvalues 
of the Hermitian matrix <I>, and A (A) = rii<j(^« ~ ^^j) ^^^^ known Vandermonde 

determinant. The /3-deformation replaces the matrix integrand by its /3 power, so the 
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partition function becomes 

Z(/3) = e^(/^) = J il[dX,)A^^iX)e-i^^^^'^\ (2.2) 

In this paper we mostly consider a cubic polynomial potential 

W{^) = \<^^ + \<^\ (2.3) 

We can compute the free energy perturbatively for small gs, by expanding the exponential 
and reduce the computation to the expectation values of multi-trace operators in Gaus- 
sian matrix model. For the correspondence with Dijkgraaf-Vafa conjecture and topological 
string theory, we consider two-cut solution of the matrix model in the large limit. The 
eigenvalues of the matrix $ distribute continuously around the two extrema of the potential 
(j2.3p . Here one of the extrema is actually a local maximum, so we need to use an analytic 
continuation to anti-Hermitian matrix, so that the expansion around the local maximum is 
the usual Gaussian matrix model. 

The perturbative calculations in the case without /3-deformation, i.e. (3 = 1, were studied 
in details in [17]. The idea is similar in the /3-deformed case, and the partition function of 
the two-cut solution can be written as a two-matrix model 

^ = N,H2.A,H2nr, I (n''«)^^''w(n''''.)^=''('')=-"'-'*-'-"''*'-"''*'*'.p-'') 

where fii and are the eigenvalues of the two Hermitian matrices $1, $2) and the potentials 
are 

W,{^,) = tr[ici>2+l(|)|$3]^ 

W2{^2) = -tr[^^i-^(|)^$i], 
= ^-/3iViiV2log(^) 

k=i P p=o 

Here the interaction term comes the exponentiation of the interaction of the eigenvalues 
fii and I'j in the Vandermonde determinant. We have rescale the matrices and <I>2 by 
a factor (^)2 to normalize the Gaussian potential, and neglect the unimportant overall 
factor in the partition function from the scaling. 

The correlators of multi-trace operators can be computed recursively by the loop equa- 
tion of the /3-deformed matrix model |22^ [23] . Denote the correlators 



/ D/3$tr(«>'=i)tr($'=2) • • • tr(«>*^™) exp(-i£|^) 
jD[s<!>eM-'-^) 



Ck,M,-,k^{N,(3) = ^ ^ J . ' ' (2-6) 
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where $ is a x Hermitian matrix, and the /3-deformed measure is Djs^ ~ ([li dXi) A'^f^ [X) . 
Here the indices ki are non-negative integers and the correlator is understood to be zero by 
convention if one of the indices is negative. For the cases of odd X^I^i the correlators 
vanish due to the reflection symmetry of the Gaussian potential. Also obvious is the situa- 
tion of one the indices ki = 0, where we can simply pull out a A'^ factor since tr(l) = A^. The 
loop equation, which is the Ward identity of matrix model, provide the following recursion 
relation for the correlators 

fcm— 2 m—1 

1=0 1=1 
+ {l-P){km-l)Ck,,..,k^.,,k^-2. (2.7) 

Here the index k^^ for the recursion should be strictly positive k^ > 0. For the undeformed 
case /3 = 1, the above recursion can be easily understood from the Feynman rule of m 
vertices, where a line of the last vertex can be contracted with another line from itself or 
a line from another vertex. From the recursion relation (j2.7p and the initial conditions 
Cq = N, one can then compute the correlators recursively |23j . 

We would like to obtain the genus expansion of the free energy. It turns out that there 
is a nice shift of the t'Hooft parameters so that the odd terms in the expansion vanish [1]. 
This is similar to the shift of mass parameters — t- + ^^^^ in Seiberg-Witten theory 
studied in p^[T3]. 

There are two contributions to the free energy, known as the non-perturbative part 
and the perturbative part. The non-perturbative contributions come from the measure of 
the matrix integral and the U{N) volume factor [26], which is obtained by evaluating the 
partition function (j2.4|) . without the interaction term VF($i,$2) and setting Qs = 0. The 
case of /3 deformation is computed in [4] . Up to an unimportant constant from the analytic 
continuation to Gaussian potential $2 — >■ ^^2, the result is 

i=l 



I-3/3 + /32 l_/3 l_5^2^^4 

where the t'Hooft parameters are defined as 

ti = gsNi, i = l,2. (2.9) 

We see naively there are odd Qs power terms in ()2.8p . In refined topological string theory, 
the free energy comes from the effective action of integrating out charged BPS particles in 
graviphoton background, and as such only even power terms appear. The situation is 
remedied by a shift 

U^k-^^g„ i = l,2. (2.10) 
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The shift cancels the odd power terms in the non-perturbative contribution ()2.8|) . We use 
the calhgraphic symbol to denote the free energy after the shift 

^n.,. = - - ^ - 5760/33 ^ + ^(5^)]- (2-11) 

Next we compute the perturbative contributions. We can neglect the unimportant terms 
in the first line in VF(^>i, $2) in ()2.5p . and keep the contributions that are perturbative for 
small Qs- We expand the remaining exponents in (I2.4p . and use the recursion ()2.7p . We find 
the first few terms 

Fj,ert = {Ni - N2)W\2Nl - + 2iV|) - 9/3(/3 - l){Ni + N2) + (5/3^ - 9/3 + 5)]|| 

+ [2/3^(8iVi^ - 91iVfiV2 + 177iVf ^1 - 91iViAf| + 8A^|) - /32(/3 - l)(59A^f - 81iVf^2 
-81iViA^| + 59Ar|) + /3(73/32 - 132/3 + 73)(A^f + _ 2/3(92/3^ - 153/3 + 92)A^iiV2 

-(/3 - l)(30/32 - 43/3 + 30)(iVi + A^2)]^ + 0{gl). (2.12) 

Again we find that the shift (j2.10p eliminates the odd Qs power terms in the perturbative 
free energy. Putting the non-perturbative and perturbative contributions (after the shift) 
together, we can extracted the higher genus refined amplitudes J^^'^'d) by the expansion 

00 

J'n.p. + T^ert = ^ g2(n+g)^2 ^l^g-2n _ i)2n_^{n,9) (^^ ^ _ (2.13) 
n,g=0 

We list the first few refined topological amplitudes. Here the unrefined cases of 
have been computed in [17], and for completeness we also display them here. 

i=l 

+177tl4 - 91titi + 84) + O(t^), (2.14) 

2 



_^(i,o) ^ -^i^^ + i^(ti-t2) + ^(97i?- 265M2 + 97ti) 
1=1 
1 



+-L(4004tf - 19401t?t2 + 19401tit| - 4004t^) + O(t^), (2.15) 
18 

^^'^'^ = -E^ + ^(ti-t2) + ^(7t?- 31^2 + 7^2) 

i=l 

+-(332tf - 2769t?t2 + 2769tit^ - 332t^) + 0{t^) (2.16) 
9 

2 

7 131 22709 1 
_^(2,o) ^ _^__ + _ + __(t^_t2) + _(581203t?- 1449550tit2 + 581203*2) 

i=l * 

23420099(tf + tj) - 100452495*1*2 (ti - ^2) .^4^ .2 17) 

120 ^ ^ ^ ' ^ 
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^ 7 17 4133, , 1 



•^^'■'^ = -Elii5^ + I^ + ^(*i-*2) + ^(120367t?- 342334iii2 + 120367*2^ 



^ 5206799(tf +ti)- 26406615tif2(ti-t2) ^ ^^^4) (2.I8) 



30 

2 



1 35 

jr(o,2) ^ _ ^ + _ + (338t? - 1632tit2 + 338ti) 



240*2 6 

+^[16533(t? + 4) - 151100tit2(ti - ^2)] + 0{t^). (2.19) 

The leading terms in from the non-perturbative contributions are basicahy the singu- 

lar terms appearing in integrating out a massless BPS particles in the general graviphoton 
background near a conifold point in refined topological string theory. We note that there 
is a (—1)" factor diff'erence comparing with the convention in |15) . probably due to some 
analytic continuation between these two kinds of calculations. 

3 Review of Dijkgraaf-Vafa geometry 

Dijkgraaf and Vafa conjectured [8] that the Hermitian matrix model with a degree n + 1 
polynomial potential VF(<I>) is dual to the topological string theory on a local Calabi-Yau 
three-fold geometry, defined by the following curve on the coordinates (n, v, y, x) 



uv 



y^ + W'{xf + f{x), (3.1) 



where f{x) is a degree n — 1 polynomial whose coefficients parametrize the complex struc- 
ture moduli of the Calabi-Yau geometry. The A-periods of the Calabi-Yau geometry are 
identified with the filling fractions, i.e. the t'Hooft parameters in large A*" limit, around the 
extrema of the potential in the n-cut solution of the matrix model. 

We specialize to the case of cubic potential W{x) in (|2.3p and introduce some notations. 
The complex deformation f{x) split the double roots (oi, 02) in the equation W'{x)'^+f{x) = 
0, whose roots are denoted as (0^,0^^,02^,0^) = (xi, X2, X3, X4). We define the complex 
parameters 

ix2 - Xi)'^ (X4 - X3)2 
^1 = ^ > ^2 - ^ , (,3.2j 

which parametrize the complex structure moduli space of the Calabi-Yau geometry. The 
moduli space and its singular divisors are studied in details in [H]. Here the divisor ziZ2 = 
behaves like the conifold divisor. There are two other singular divisors 

I{zi,Z2) = ^[(X3 + X4) - (Xi + X2)]^ = 1 - 2(zi + 2:2), 
J(zi,22) = (Xi - X3)(X2 - X3)(xi - X4)(X2 - X4) 

= 1 - 6{zi + Z2) + 9zl + UziZ2 + 9zl (3.3) 

where the singular divisor J = also behave like the conifold divisor, and the higher genus 
topological string amplitudes expended around the divisor t will have gap like behavior [HI 
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116] . On the other hand, the singular divisor I = is actually an essential singularity where 
the perturbative solutions to the Picard-Fuchs equation have zero radius of convergence. 

The planar free energy J^^^^ of the matrix model are determined by the well known 
equation ^^-^ = Hj, where the periods are 



1 . , „ 1 



A 



as contour integrals over cycles of x-plane, with the one-form differential Xdx = \JW'{x)'^ + f{x)dx. 
The expansion for the A-periods to the first few orders are 

ti = J - |-(2^i + 3z2) - |^(4z? + I3ziz2 + 9zl) + 0{z*), 

t2 = -^ + ^i2z2 + 3zi) + ^i4zi + l3ziZ2 + 9zl) + Oiz^). (3.5) 

To compute the period integrals, one finds differential operators C{zi, Z2) such that CX 
is a total derivative of x variable. One also needs to be careful about possible residua, 
because of which an operator C may not annihilate the periods even though CX is a total 
derivative [14] . More precisely, up to the second order differentials, there are three linearly 
independent operators whose actions on A are total derivatives, which are the followings 



C2 = 



= [(3 - 2zi - 6z2)d,, - 2zi(l - 2zi - 6z2)dl + (zi o Z2)] 

+2(1 - 5zi - 5z2 + 6zf + AziZ2 + 6z^)d^^d^^, 
- -6 + 36(zi + Z2) - Q{9zl + UziZ2 + 9z1) + 2{zi + Z2)[-l + Qzi + Qz2 - llzf - lOzi^s 
+ 6(zi - Z2fizi + Z2)]d,,d,^ + [{7zi - 3z2 - 39zf - lSziZ2 + 9z| + 46z? 
+mzlz2 + 26ziz| - 6z|)a^i + 2zi(l - Qzi - 6z2 + V3z\ + \\zxZ2 + 5z| - lO^f 

-^Z\Z2 + lOzizf + ^zl)dl^ + (Zi O Z2% 

[(5zi + 3z2)d,^ + 2zi(2 - 5zi - 3Z2)5,', - (^i ^ ^2)] 

-2(zi - Z2)(l - 3zi - 3z2)a,,a,, (3.6) 



All three operators annihilate the A-periods tj. After adding proper constants from regu- 
larizing the integrals, the B-periods Hj can be annihilated by the operators L\ and £2, but 
not by £3. 

The three point couplings are rational functions of z\^Z2, and have been computed in 



_ 1 - (621 + 5^2) + 3(3z\ + 3zxZ2 + 2zl) 



\^I{Z\,Z2) 

1 — 321 — 5Z2 
16I(zi,Z2) 



and the others are related by symmetry. 



6 



4 Genus one formulae 



For the genus one case, we have the following exact formulae from refined topological string 
^(0,1) = _liog(det(|^))-llog(ziZ2)-^logI(zi,^2) + ^logJ(^i,^2), 

= -^l0g{z,Z2)-^l0gJizi,Z2), (4.1) 

where /, J are the singular divisors in ()3.3p . These formulae can be checked perturbatively 
with the expansions ()2.15p . ()2.16p . For the case of amplitude, the exact formula can 

be also derived from the loop equation in matrix model [5]. However, the case of T^^'^^ 
amplitude, which comes from the /3-deformation, seems much more challenging to obtained 
exactly from the matrix model methods. 

In [1], the perturbative expansion for J^(^'^^ amplitude is used to test the idea of quantum 
integrability in the Nekrasov-Shatashvili (NS) limit [251 I21j . Here we should follow the 
method in [12], and check the formula exactly. First we define the deformed periods 



ti= f \dx, Hj = * Xdx, (4.2) 
JAi JBi 

where the deformed one-form differential Xdx is determined by the quantizing the curve 

pH{x) = [W'{xf + f{x)]^{x), (4.3) 

where ^{x) = exp(i Xdx) is a wave function. The canonical quantization relation is 
imposed so that the operator p = edx- We can use the WKB expansion of the wave 
function 

oo 

A(x) = ^A„(x)e-, (4.4) 
and solve perturbatively for the first few terms as the followings 



A(x) = Ao(x) = ^JW'{xY + f{x), 

X'{x) , . X _ 2A"(x)A(x) - 3A'(x)2 
' X{x) ' ^^^''> ~ 8A(x) 



Ai(x) = -2-;7-^, A2(x) = ^,^^.3 , ••• (4.5) 



The odd terms are total derivatives and vanish when integrating over a contour. We are 
interested in the even term contributions. To compute the contour integrals, we should 
relate the integrand X2n{x) to the action of some deferential operators on A(x) plus some 
total derivatives of x. Such a differential operator is found in [l], whose derivatives are with 
respect to the roots Xi{i = 1,2,3,4) of the quartic equation W'{x)'^ + f{x) = 0. Here we 
should write the operator as derivatives of complex structure parameters zi,Z2, which is 
more convenient for calculations. Furthermore, we find that because of the extra Picard- 
Fuchs operator with non-vanishing residue, the operator is not uniquely determined but is 
ambiguous by addition of an operator whose action on A(x) is a total derivative. 



7 



We consider the differential operator without the constant and the dz^dz2 terms, which 
can be always eliminated using the Picard-Fuchs operators in (13. 6p . Up to a total derivative 
denoted as ^(x), we find A2 = T>2X + h'2{x), where the operator 

V2 = l/[6J(zi, Z2)^(l - 5zi - 5z2 + 6zf + AziZ2 + 6z|)]{[l + 12(zi - Z2) - 2(121zf + sOziZa 
-2f>zl) + 4(33l2? + 543zf Z2 + 229ziz| - 15zf ) - (3303^f + 8084zf Z2 + 7978z?z| 
+2772ziz| + 1354) + 16(243zf + 756z^Z2 + 1123z?4 + 909z?4 + 270zi4 + 27^) 
-12(144z? + 5252^22 + 967z^z| + 1314zf z| + 894z?z| + 225zi4 + 27z^)]a^i 
+ [-1 + 2(13zi + 7^2) - 2(100z? + 141ziZ2 + 39^) + 4(169z? + 388z?Z2 + 285zi4 
+54^) - 3(3534 + 1098zf Z2 + 1356z?z| + 678ziz| + 99^) + 6(105^ + 403z^Z2 
+694zf 4 + 59444 + 225zi4 + '^^4)]dl^ } + (^1 ^ Z2). (4.6) 

There is still an ambiguity due to linear combination of the Picard-Fuchs operators C2 and 
C^. We will consider the NS limit and use the perturbative expansion of J^^^''^) amplitude 
()2.15p , in order to fix the ambiguity. The operator P2 turns out to be the correct operator 
which compute the deformed period from the leading term. So we find 

ti = ti + {V2ti)e^ + 0{e^), 

n, = ni + iV2U,y + 0{e^). (4.7) 

We denote the deformed prepotential in the NS limit as 

00 

^W(t,) = E-^("'°)(ti)(J)", (4.8) 

n=0 

where we use a normalization factor of 4" for convenience. As in |12| we expand the equation 
for J^^^^ in terms of deformed periods, to derive differential equations for higher genus terms 
in the NS limit 

= a,-/W(t,)-n, 

= 9^,^(0,0) _ n. + ,2^du^f^ ^ (P2t,)(5t,9i^..F(0'0)) - V2U,] + 0{e^). (4.9) 

After some algebra, the order equation can be shown to be equivalent to 

= (4.10) 

where Cz^ZjZk the three-point Yukawa couplings, and Pzjz^ the coefficients dz^dz,. in 
the differential operator P2- We check the above equation (|4.1U|) is satisfied exactly, using 
the three point functions (j3.7p and the differential operator D2 in (|4.6p . 

5 Higher genus 

The refined topological string amplitudes at genus n + g >2 can be written as polynomials 
of the propagators 5^'^^ of degrees 2n + 3^ — 3, with rational functions of complex structure 
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moduli zi^2 as coefficients. Here for the local geometry we can set the Kahler potential to 
be constant and only need the double-index S'^'^j propagators. The propagators S^^^^ are 
defined in [3] by their anti-holomorphic derivative relation with the three point functions, 
which can be integrated to 

^^'^^^w. = -nu+/S.p (5.1) 

where the fzlzi are holomorphic ambiguities from the integration, rational functions of zi^2, 
and have been determined in [14j . The holomorphic derivatives of the propagators form a 
closed algebra 

dz,S'^'^ = Cz,z^z„S'^'-S'^'" - /S.^^^^"™ - fziz^S'-'- + hl'\ (5.2) 

where hl]^^ are also rational functions of zi^2- We also fix them and find it is sufficient to 
use J{zi,Z2)^ as the denominator. 

We integrate the refined holomorphic anomaly equation and fix the holomorphic am- 
biguity by boundary conditions near the singular divisors of the Dijkgraaf-Vafa geometry. 
Here the ansatz for holomorphic ambiguity at genus g is 

f _ K,g{Zl,Z2) 

where hn^g{zi, Z2) is a symmetric polynomial of zi and Z2, and it turns out that for low 
genus, its degree is lower than 9(n + g — 1) [16] . 

To utilize the gap condition near the conifold divisor J, we choose a point on the divisor 
from its intersection with the line zi = Z2- Near the point (21,^2) = (|) |)) the good local 
coordinate is 

Zc,l = Zi- Z2, Zc,2 = 1 - 4(21 + Z2). (5.4) 

There are 3 power series solutions to the Picard-Fuchs equation near this point 

Wl = Zc^i^l + Zc,2, 

W3 = (44l+3Zc,l^c,2) + (184l^C:2 + ^^c,l42) + 0(^c)- (5-5) 

One choose the first two solutions as the flat coordinates tc,i = Wi near this point [16j. 

The singular terms of the refined topological string amplitudes expanded near the divi- 
sors 2:12:2 = and J = are 

*2 

•^^"'^^ = ^ ^2,-2 +Oitl„tl2)- (5.6) 
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Here the factor of (— 2)^^^^"^) is due to a normalization of the flat coordinate tc,2i and the 
constants c"'^ appear e.g. in (12.171 12. 18^ I2.19P for genus two. These gap conditions (15. 6p 
are sufficient to fix the holomorphic ambiguities, since any non-zero ansatz in ()5.3p with 
hn,gizi, Z2) a polynomial of degree less than 12(n -\- g — 1) can not cancel all the poles in 
the denominators, so would affect the singular terms. 

We find the exact refined formulae for the genus 2 amplitudes and check the expansions 
near the point {zi, Z2) = (0, 0) agree with the higher order terms from perturbative matrix 
model calculations in (j2.17| [27T81 12.19p . For example, the formula for is 

= 1 1 .0 ^2 2 ,2 ["i^'^'^ + 2QiQ2g^^^^ + alS^-^-] + /2.0, (5.7) 
Here the coefficients of the propagators are 

oi = 1 - ISzi - 6z2 + 452;f + 422122 + Qzf, 02 = 01(21022), (5.8) 
and the holomorphic ambiguity is 

/2,o = J^2j4 [-(11^1 + ^2 + ll^i) + (231^? + 3892^22 + •••)- (20792^ + 55472322 

+mAzlzl + •••) + 15(6932^ + 24872^22 + 23242^21 + •••)- 9(34652^ + 154052^22 
+194352^2^ + 83822^2^ + •••) + 3(187112^ + 980372^22 + 1727912^2^ + 217572^2^ + • • • ) 
-3(187112^ + 1116992^22 + 2662742^ 2^ + 509732^2^ - 2727222^2^ + • • • ) 
+27(21 - 22)^(8912^ + 76952^22 + 326252^ 2I + 673332^2^ + •••)], (5-9) 

where the • • • denote the terms implied by the symmetry 21 o 22. 
Acknowledgments 

We thanks Albrecht Klemm for discussions and collaborations on related papers. MH is 
supported by the "Young Thousand People" plan by the Central Organization Department 
in China. 

References 

[1] M. Aganagic, M. C. N. Cheng, R. Dijkgraaf, D. Krefl and C. Vafa, "Quantum Geometry 
of Refined Topological Strings," JHEP 1211, 019 (2012) [arXiv: 1105.0630 [hep-th]]. 

[2] G. Akemann, "Higher genus correlators for the Hermitian matrix model with multiple 
cuts," Nucl. Phys. B 482, 403 (1996) [hep-th/9606004]. 

[3] M. Bershadsky, S. Cecotti, H. Ooguri and C. Vafa, "Kodaira-Spencer theory of gravity 
and exact results for quantum string amplitudes," Commun. Math. Phys. 165, 311 
(1994) [hep-th/9309140]. 

[4] A. Brini, M. Marino and S. Stevan, "The Uses of the refined matrix model recursion," 
J. Math. Phys. 52, 052305 (2011) [arXiv:1010.1210 [hep-th]]. 



10 



[5] L. Chekhov, "Genus one correlation to multicut matrix model solutions," Theor. Math. 
Phys. 141, 1640 (2004) [Teor. Mat. Fiz. 141, 358 (2004)] [arXiv:hep-th/0401089]. 

[6] L. Chekhov and B. Eynard, "Matrix eigenvalue model: Feynman graph technique for 
all genera," JHEP 0612, 026 (2006) [math-ph/0604014]. 

[7] L. O. Chekhov, B. Eynard and O. Marchal, "Topological expansion of /3-ensemble 
model and quantum algebraic geometry in the sectorwise approach," Theor. Math. 
Phys. 166, 141 (2011) [arXiv: 1009.6007 [math-ph]]. 

[8] R. Dijkgraaf and C. Vafa, "Matrix models, topological strings, and supersymmctric 
gauge theories," Nucl. Phys. B 644, 3 (2002) [hcp-th/0206255]. 

[9] B. Eynard, M. Marino and N. Orantin, "Holomorphic anomaly and matrix models," 
JHEP 0706, 058 (2007) [hep-th/0702110 [HEP-TH]]. 

[10] B. Eynard and N. Orantin, "Invariants of algebraic curves and topological expansion," 
Commun. Num. Theor. Phys. 1, 347 (2007) [math-ph/0702045]. 

[11] B. Haghighat, A. Klemm and M. Ranch, "Integrability of the holomorphic anomaly 
equations," JHEP 0810, 097 (2008) [arXiv:0809.1674 [hep-th]]. 

[12] M. -X. Huang, "On Gauge Theory and Topological String in Nekrasov-Shatashvili 
Limit," JHEP 1206, 152 (2012) [arXiv:1205.3652 [hep-th]]. 

[13] M. -X. Huang, A. -K. Kashani-Poor and A. Klemm, "The Q deformed B-model for rigid 
^f = 2 theories," Annales Henri Poincare 14, 425 (2013) [arXiv:1109.5728 [hep-th]]. 

[14] M. X. Huang and A. Klemm, "Holomorphic anomaly in gauge theories and matrix 
models," JHEP 0709, 054 (2007) [arXiv:hep-th/0605195]. 

[15] M. -X. Huang and A. Klemm, "Direct integration for general Q backgrounds," Adv. 
Theor. Math. Phys. 16, no. 3, 805 (2012) [arXiv:1009.1126 [hep-th]]. 

[16] A. Klemm, M. Marino and M. Ranch, "Direct Integration and Non-Perturbative Effects 
in Matrix Models," JHEP 1010, 004 (2010) [arXiv:1002.3846 [hep-th]]. 

[17] A. Klemm, M. Marino and S. Thciscn, "Gravitational corrections in supersymmetric 
gauge theory and matrix models," JHEP 0303, 051 (2003) [arXiv:hep-th/0211216]. 

[18] D. Krefl and J. Walcher, "Extended Holomorphic Anomaly in Gauge Theory," Lett. 
Math. Phys. 95, 67 (2011) [arXiv: 1007.0263 [hep-th]]. 

[19] D. Krefl and J. Walcher, "Shift versus Extension in Refined Partition Functions," 
arXiv: 1010.2635 [hep-th]. 

[20] D. Krefl and J. Walcher, "ABCD of Beta Ensembles and Topological Strings," JHEP 
1211, 111 (2012) [arXiv: 1207. 1438 [hep-th]]. 



11 



[21] A. Mironov and A. Morozov, "Nekrasov Functions and Exact Bohr-Sommerfeld Inte- 
grals," JHEP 1004, 040 (2010) [arXiv:0910.5670 [hep-th]]. 

[22] A. Mironov, A. Morozov, A. Popolitov and S. .Shakirov, "Resolvents and Seiberg- 
Witten representation for Gaussian beta-ensemble," Theor. Math. Phys. 171, 505 
(2012) [arXiv: 1103.5470 [hep-th]]. 

[23] A. Morozov and S. Shakirov, "The matrix model version of AGT conjecture and CIV- 
DV prepotential," JHEP 1008, 066 (2010) [arXiv: 1004.2917 [hep-th]]. 

[24] N. Nekrasov and A. Okounkov, "Seiberg-Witten theory and random partitions," 
arXiv:hep-th/0306238. 

[25] N. A. Nekrasov and S. L. Shatashvili, "Quantization of Integrable Systems and Four 
Dimensional Gauge Theories," arXiv:0908.4052 [hep-th]. 

[26] H. Ooguri and C. Vafa, "World sheet derivation of a large N duality," Nucl. Phys. B 
641, 3 (2002) [hep-th/0205297]. 



12 



